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Abstract. We provide a framework to build periodic boundary conditions on 
the pseudosphere (or hyperbolic plane), the infinite two-dimensional Riemannian 
space of constant negative curvature. Starting from the common case of 
periodic boundary conditions in the Euclidean plane, we introduce all the needed 
mathematical notions and sketch a classification of periodic boundary conditions 
on the hyperbolic plane. We stress the possible applications in statistical 
mechanics for studying the bulk behavior of physical systems and we illustrate how 
to implement such periodic boundary conditions in two examples, the dynamics 
of particles on the pseudosphere and the study of classical spins on hyperbolic 
lattices. 



1. Introduction 

The use of periodic boundary conditions has become a standard tool in statistical 
physics, and more particularly in computer simulations, to extract from the behavior 
of finite systems relevant information on the properties in the thermodynamic limit. 
In Molecular Dynamics and Monte Carlo simulations, the number of particles, or 
more generally degrees of freedom, that can be studied with present-time computers 
remains limited (say, of the order of 10 6 or less) and extrapolating to the "bulk" 
behavior of a macroscopic system requires to somewhat minimize the influence of 
the system's boundaries. This is commonly done by employing periodic boundary 
conditions [UEJE]. 

The practical implementation of such boundary conditions is well developed in 
Euclidean, "flat" spaces. However, much less has been achieved with regard to curved 
spaces. We do not consider here spherical-like spaces characterized by a positive 
Gaussian curvature, such as the surface of a sphere in three-dimensional Euclidean 
space, because such spaces have a finite extent; as a result, the thermodynamic limit 
is only attained by letting the curvature go to zero. We rather focus on hyperbolic 
spaces, characterized by a negative curvature, which can indeed be infinite. Due to 
the peculiar character of the hyperbolic metric, the thermodynamic limit in the case 
of an open space may crucially depend on the nature of the boundary conditions: 
for a homogeneous, simply-connected hyperbolic space, the perimeter (or area in 3 
dimensions, etc) of the boundary of a disk of radius r is proportional to the area of 
the disk (volume in 3D, etc), when r is much larger than the radius of curvature. 

Negatively curved spaces show up in different areas of physics. The more 
prominent one is general relativity and cosmology where the curvature of spacetime 
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is a fundamental ingredient. For the same reason, and usually motivated by quantum 
gravity and string theories, quantum field theories have been studied in curved 
spacetime [4 , in particular spacetimes with a hyperbolic spatial section [U [5j |6j [71 El E] • 
Negatively curved spaces have also recently been considered in the context of statistical 
physics and condensed matter theory. There, the motivations are both practical 
(describing the behavior of newly designed mesoscopic and nanoscopic objects with 
exotic shapes, curvatures, and topological properties p~Ql \TT\ \T2\ [321 HH E]) and 
theoretical (understanding how curvature influences the critical behavior and the phase 
transitions of classical statistical systems [5| t8| fTB l fTT l fT8| fT9l f2Q | f2T | f22| 123 ^ f2^ 1 f25| 
[26l [271 EH [29] ) . Still in condensed-matter theory, negatively curved spaces appear in 
studies of the Quantum Hall effect [TQl [TTJ [12] and in the framework of "geometrical 
frustration" [30l [311 [32]. In the latter case, the hyperbolic reference space serves 
either for providing crystal-like templates to amorphous solids [33] or exotic liquid 
crystalline structures (such as the "infinite periodic minimal surfaces" associated with 
the bicontinuous cubic phases of liquid crystals formed by amphiphilic molecules and 
water [34]), or for building toy models of supercooled liquids and metallic glasses 
[T6l [3TJ [35] . Finally, the free motion of mass points on hyperbolic surfaces has been 
used as a prototype to study classical and quantum chaos [36], and properties of 
random walks have also been investigated [371 EH] • 

In the above listed studies, need for periodic boundary conditions may arise for 
different reasons, (i) First, and as mentioned above, they may be useful to simulate, 
as closely as possible, the bulk behavior of a macroscopic sample through the study 
of finite-size models, (ii) Secondly, periodic boundary conditions, as will be further 
detailed in the following, are intimately connected with tessellations of space, i.e., 
tilings of space by identical replicas of a unit cell; as a consequence, their study has a 
direct bearing on the characterization of crystalline structures, and thus, in the present 
case, on "hyperbolic crystallography", (iii) From a more topological point of view, 
periodic boundary conditions are also related to the generation of compact manifolds 
with nontrivial topology. Such multi-connected hyperbolic manifolds are for instance 
relevant to "cosmic topology", which investigates the possibility that the universe is 
not simply connected [39 , 4(1, and to studies on classical and quantum chaos [36] . 

In the following, we restrict ourselves to the consideration of periodic boundary 
conditions on the "hyperbolic plane" or "pseudosphere" (also called "Bolyai- 
Lobachevsky plane"), which is the infinite, simply-connected two-dimensional space 
with constant negative curvature [4TJ |42]. Our goal is to provide a framework for 
implementing periodic boundary conditions with chosen properties. Such properties 
will generally follow from physical requirements and may be specific to the problem 
under study. 

The rest of the article is organized as follows. In section II we revisit the case of 
periodic boundary conditions in the Euclidean plane in order to introduce the two- 
fold approach, geometrical and topological, to such conditions. We provide in the next 
section the basic mathematical formalism for studying periodic boundary conditions 
in the hyperbolic plane. In section III, we address the classification of the possible 
periodic boundary conditions in the hyperbolic plane and give a guideline for building 
the proper conditions adapted to a given physical problem. Finally, in section IV, 
we illustrate how to implement the proposed formalism by considering two examples 
of statistical mechanical models: the dynamics of particles on the pseudosphere and 
spin models on hyperbolic lattices. We complement the presentation by appendices 
in which we summarize basic elements on hyperbolic geometry (Appendix A), on 
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Figure 1. Square periodic boundary conditions in the Euclidean plane, (a) 
Square with pairing of the facing edges, (b) One-hole torus obtained by gluing 
paired edges of the square. The embedded graph (black and dotted lines) 
corresponds to the paired edges. 



tessellations (Appendix C). 



isometries and Fuchsian groups (Appendix B), and on hyperbolic trigonometry and 



2. The Euclidean case revisited 

We first consider the well-known case of the Euclidean plane E 2 . Implementing 
periodic boundary conditions in E 2 consists in choosing a primitive cell (which contains 
the physical system) such that it can be infinitely replicated to tile the whole plane 
with no overlaps between copies of the cell and no voids. To ensure smoothness and 
consistency, the edges of the primitive cell have to be paired in a specific way: leaving 
the cell through one edge implies to come back by another edge, which should be 
equivalent to exploring the whole tiling of the plane. 

The simplest periodic boundary conditions in E 2 are the "square" ones: taking 
a square as the primitive cell, one can obviously replicate it to form a square lattice, 
which corresponds to a particular tiling of the plane, and one can pair facing sides 
of this primitive square (see figure [T(a)) ) . This operation can be framed in a more 
mathematical setting [43]. Indeed, in the language of group theory, periodic boundary 
conditions are just an equivalence relation. In the present case, choosing any two 
points in E 2 , i.e., any two pairs {x\,yi) and (#2,2/2) m one can introduce the 
following equivalence relation: 



0i,2/i) ~ (>2,2/2) 



if 



x 2 



x\ + a n x 



V2 Vi 



(i) 



G Z and a the side length of the square. The quotient space R 2 / ~ 
associated with this equivalence relation (in geometric context, "orbit space" is 
sometimes used in pl ace of "quotient space") is then the one- hole torus, noted T 2 . As 
can be seen in figure 1(b) , T 2 can be represented by the primitive cell with identified 
edges (which is then called in two dimensions the fundamental polygon): by "gluing" 
paired edges together, one can indeed obtain a torus which is t opolo gical^ equivalent 
to the one-hole torus T 2 . Note that the torus shown in figure 1(b) is represented as 
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the projection of an object that seems to be physically embedded in the 3D Euclidean 
space E 3 and appears, with the natural Euclidean metric, curved like a doughnut; this 
is not the case of the "flat" two-dimensional torus T 2 , whose actual visualization is 
then more elusive. 

The way to pair the edges is crucial. Other pairings of the edges of the primitive 
square are in fact possible. They lead to different 2D manifolds, which can be 
considered as more "exotic" than the torus, being nonorientable and intersecting with 
themselves in E 3 , such as the Klein bottle and the projective plane (this latter, in 
addition, having a positive curvature) [42 , 43 . Actually, the way to pair edges reflects 
the equivalence relation underlying the periodic boundary conditions. 

As in the simple example above, we will use in this article, depending on the 
context, the two views concerning periodic boundary conditions: 

(1) the "geometrical" one, based on the tiling of the infinite and simply connected 
space X endowed with its metric (E 2 in the above example) by identical 
(congruent) replicas of the fundamental cell; such a tiling or tessellation involves 
a discrete subgroup T of isometries of the space X (i.e., displacements preserving 
the metric), whose elements transform the replicas of the cell one into another. 

(2) The "topological" one, using multi-connected manifolds (here, surfaces) to 
represent the quotient space, X/Y. The link between Y and the equivalence 
relation ~ defined in equation ([!]) is that two points x, y of X are equivalent if 
and only if there exists an element 7 of the group Y such that j(x) = y. For 
example, in the case of equation [I] Y is generated by the two translations T a £ and 

-Lay 

These two points of view are fully compatible and closely connected. A "visual" link 
can be made between them by gluing paired edges of the fundamental polygon (here, 
the square with identified edges) to obtain a manifold representing the quotient space 
(here, the one- hole torus), with the glued edges forming a graph embedded in the 



manifold as shown in figure 1(b) As we shall see below, this latter aspect is essential 
for more general geometries. 

In the 2D flat space, alternatives exist to build periodic boundary conditions. 
One can consider a hexagonal cell, also with paired opposite edges (see figure [2(a)] ) • 
In this case, the underlying equivalence relation for two pairs (x\,yi) and (#2? 2/2) m 
R 2 is: 

3a 

a x 

(2) 



(2:1,2/1) ~ (2:2,2/2) if < 



X2 = x 1 + — n x 
2/2 = 2/1 + ~^r- % 



where n x , n y and a are defined as below ([!]). The associated quotient space is also the 
one- hole torus T 2 , but the graph formed by glued edges of the hexagon is emb edded 
in T 2 in a totally different way than that used before, as illustrated in figure 2(b) 



As in the case of the square, other pairings of the edges are possible. However, the 
resulting manifolds associated with the quotient space are again exotic, nonorientable 
surfaces. 

In what follows, we generalize the present considerations to the hyperbolic plane. 
As they are simpler and usually more physical, we only consider periodic boundary 
conditions that can be built from a primitive cell which is a regular polygon with a 
finite area (due to the pairing requirement, the number of edges must in addition be 




Figure 2. Hexagonal periodic boundary conditions in the Euclidean plane, (a) 
Hexagon with pairing of the facing edges, (b) One-hole torus obtained by gluing 
paired edges of the hexagon. The embedded graph (black a nd do tted lines) 
corresponds to the paired edges. Note the difference with figure [T(b)] 



even). The associated regular tessellation of the plane is generically noted {p, q}, where 
p is the (even) number of polygon edges and q the number of polygons that meet at a 
vertex. (We shall discuss the constraints on p and q below.) We also avoid "twisted" 
boundary conditions and therefore restrict pairings of the sides of the primitive cell 
to schemes that lead to quotient spaces representable as orient able surfaces. With the 
above restrictions, the surface is then closed and of finite extent with respect to the 
metric under consideration, namely, compact. 



3. Generalization to the hyperbolic plane: fundamentals 

3.1. Tessellations, quotient spaces, and fundamental polygons 

Going back for a moment to the case of the Euclidean plane, we recall that well 
known geometric constraints limit the possible regular tilings {p, q] to those verifying 
(p— 2)(q— 2) = 4. This includes the tilings we have already considered, {4, 4} (tiling by 
squares) and {6, 3} (tiling by hexagons), plus a third one, the tiling {3, 6} by triangles 
which is dual to the {6,3}. This {3,6} tiling however involves a regular polygon 
with an odd number of sides, which cannot be paired to generate periodic boundary 
conditions. (It is nonetheless possible to build periodic boundary conditions by using 
the {3, 6} tiling, but in that case the fundamental polygon is made of two triangles 
joined to form a parallelogram with paired opposite sides; but as we have stressed 
above, we only consider in this article cases where there exists a regular fundamental 
polygon.) 

In the hyperbolic plane, which we shall denote H 2 irrespective of the model used 



for its representation (see below and Appendix A), regular tilings exist provided 



(p-2)(g-2)>4, (3) 

which now leads to an infinite number of possibilities [41] |42] |44] . The richness of the 
tessellations in hyperbolic space, compared to that in Euclidean space, can be traced 
back to the existence of an intrinsic lengthscale, the radius of curvature ft -1 , where 
— ft 2 is the Gaussian curvature of H 2 [36] . 
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Another property brought about by hyperbolic geometry is that quotient spaces 
associated with periodic boundary conditions have a greater variety than in Euclidean 
space. In the latter, we have seen that E 2 /Y is a one-hole torus, which is a closed 
orientable surface of genus 1. (The genus of a connected, orientable surface is the 
number of holes, or handles, on it: for instance, a sphere and a disk have genus zero, 
whereas a one- hole torus has genus one.) Using the Gauss-Bonnet theorem which links 
the area A of a surface to its genus g, one finds, for any surface embedded in H 2 with 
negative curvature — ft 2 , 

A = 4iTK- 2 (g-l), (4) 

which implies that any such surface has a genus g ^ 2. As a result a (finite-sized) 
fundamental polygon in H 2 must have a genus g ^ 2 and so does the associated 
(compact) quotient space. As we shall see, there is no other restriction on the value 
of the genus in this case. 

To proceed further in our investigation of periodic boundary conditions in the 
hyperbolic plane, we focus on the way to obtain fundamental polygons. These 
polygons are direct ly connected to the discrete subgroups T of the group of isometries 
of H 2 (see [|] and Appendix B). There are a number of restrictions on T. First, 



being interested in quotient spaces H 2 /Y that are orientable manifolds, we only 
consider groups containing orientation-preserving transformations. Such groups Y 



are called Fuchsian groups [45J [46] (see Appendix B). Fuchsian groups play for the 



hyperbolic geometry a role similar to that of crystallographic groups for the Euclidean 
geometry. Note however that a Fuchsian group is non-Abelian. (Discrete subgroups 
containing elements which do not preserve orientation also exist and are generically 
called non-Euclidean crystallographic groups.) Actually, we will be interested only by 
"purely hyperbolic" Fuchsian groups which only contain transformations acting on the 
hyperbolic plane without fixed points and which do not lead to infinite polygons, i.e., 
generalizations of the Euclidean translations used in standard crystallography [45) [46]. 

To each Fuchsian group Y is associated one or more fundamental polygons. In 
general indeed, the fundamental polygon is not unique. Among the possible ones, 
two types are worth noticing: the standard or canonical fundamental polygon and 
the metric fundamental polygon. The former is defined by a side pairing which is 
given by a specific sequence of the generators of the associated Fuchsian group [§] 
The latter is defined as the Voronoi cell (in physicists' language) or Dirichlet polygon 
(in mathematicians' language) built around a suitable base point and endowed with a 
proper side pairing; more precisely, the Voronoi construction is applied around the base 
point by considering all the images of the point obtained by applying the generators 
of the relevant Fuchsian group. It has been proven that for each T, there exists a 
unique standard polygon [47] [48], but a metric one may not necessarily exist. (On 
the other hand, once the metric fundamental polygon is known, one can construct the 
associated standard polygon [45].) Note also that these two polygons usually have a 
different number of sides. 

The difference between the two types of fundamental polygons can be illustrated 
by taking again the example of the Euclidean plane. In the case of the {6, 3} tiling, 
the metric fundamental polygon is a hexagon (the corresponding Voronoi polygon is 

X Two points are images one of each other by a transformation of a subgroup T if and only if they 
satisfy the equivalence relation associated with the quotient space H 2 /V. 

§ The specific sequence of generators can be written 71 72 7f 1 T^ 1 • • • l2g-i 12 9 l^g-l ^2g = ^> 
where g is the genus and the 7-i's and their inverses are the generators of the Fuchsian group. 
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built around the center of the hexagon and not, as often in crystallography, around 
one of the vertices), whereas the standard fundamental polygon is a parallelogram. 
In the latter case, the pairing of the sides and the associated graph are similar to 
those of the square periodic boundary conditions shown in figure [I] The equivalence 
between the two fundamental polygons can be seen in terms of homotopy, as curves 
formed by paired edges of the hexagon (see figure [2(b)] ) and embedded in the torus are 
homotopic to the ones formed by the parallelogram (see figure 1(b)). The two nodes 
of the graph of figure 2(b) can indeed be assembled in one node to give the graph 
of figure [T(b)| If one now considers the {3,6} tiling, dual to the {6,3}, the metric 
fundamental polygon does not exist, as the unit cell is a triangle and has an odd 
number of edges. However, as seen above, the standard fundamental polygon exists 
and is a parallelogram. Finally, for the {4, 4} tiling, the metric and the standard 
polygon are identical and both realized by a square with paired sides. 

The convenience of studying metric polygons is that the link between the polygon 
and the associated tiling is clear from a geometrical point of view, whereas standard 
fundamental polygons are not trivially related to their corresponding tiling. Moreover, 
for a regular tiling, the associated Voronoi polygon is always regular, whereas standard 
polygons are regular only in some specific cases. In any case, as we are interested by 
periodic boundary conditions built from a primitive cell which is a regular polygon, 
we only consider metric fundamental polygons in the following. 



3.2. Metric fundamental polygon and associated graph 

Like any fundamental polygon, a metric one is the unit cell of a tessellation of H 2 by the 
action of the associated Fuchsian group V. This unit cell is usually called "fundamental 
domain", but we shall not use this terminology to avoid confusion with the concept 
of "fundamental polygon" which, contrary to that of "fundamental domain" , implies 
a specific pairing of the edges. Being a Dirichlet domain, the metric fundamental 
polygon is convex, ie, the geodesic joining any two points of the polygon is in the 
interior of the polygon and the edges are geodesic arcs. The edges are paired together, 
which means that each edge of the polygon can be transformed in the associated one 
by the action of some element 7 of the group V. As a result, the number of edges 
must be even. These transformations 7 actually form a set of generators of the group 
T [44j |45[ 06] . Since we restrict ourselves to orientation-preserving groups that act 
without fixed points, the generators are (hyperbolic) translations. If 27V is the number 
of edges of the metric fundamental polygon, there are 2N translations 7^ by which 
edge i is transformed into its associated counterpart. (Note that the labeling of the 
sides can be taken either clockwise or counterclockwise.) These generators must then 
satisfy the constraint 

7172 . . -7i . . -72AT = I- 

Note that for each 7^, its inverse 7" 1 is also present in the 2N generators. In the 
hyperbolic plane, these generators do not commute, and their associated Fuchsian 
group is non-Abelian. Concrete examples will be given below. 

In an attempt to classify the possible periodic boundary conditions, one may first 
consider the properties of the primitive cell (which we take as a regular polygon) in 
which the physical system will be placed: more specifically, its surface area and the 
number of its sides. From the Gauss-Bonnet theorem, see equation Q, the area is 
fixed by the curvature of the embedding hyperbolic space H 2 and by the genus g of the 
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quotient space H 2 /Y associated with the fundamental polygon. Notice that, the genus 
being an integer, the accessible values of the fundamental polygon areas form a discrete 
set. To make progress in building a classification, it is now convenient to focus on the 
graph formed by the identifying (or in more pictorial terms, gluing together) the edges 
of the fundamental polygon. This graph, which is embedded in the quotient space, 
contains all the information about the fundamental polygon, and a closed walk on the 
graph is equivalent to going around the boundary of the corresponding fundamental 
polygon. The constraints on such a closed walk will be detailed in the next subsection. 

The graphs can be first classified by their "genus" . The genus of a graph is defined 
as the minimum integer g such that there exists a closed orientable surface of genus 
equal to g in which the graph can be embedded without crossing itself. For instance, 
the graphs of figures 1(b) and 2(b) have genus 1 as they live on the one- hole torus 
without crossing themselves; on the other hand, a standard planar graph has genus 
as it can be embedded in the sphere without crossing itself. The genus of a graph and 
that of the associated fundamental polygon are equal by construction. 

The graphs are also characterized by the number of their vertices v, the number 
of their edges e, and that of their faces /. Euler's formula relates v, e, / and the genus 
g of the graph: 

v-e + f = -2(g-l). (5) 

Faces of a graph are regions bounded by edges. In the present case, there is one 
and only one face, corresponding to the interior of the polygon (which is a simply- 
connected surface). The quotient space is indeed entirely covered by the surface of 
the polygon whose edges form the graph. Stated otherwise, the compact manifold 
formed by the quotient space leads back to the fundamental polygon when cut along 
the edges of the graph. 

Inserting this property, / = 1, in the above equation gives the following relation 
between e and v: 

e = v + 2g-l. (6) 

The smallest number of vertices in the graph is v = 1, which implies that the smallest 
number of edges is e = 2g. 

To further investigate the constraints set on the characteristics of an acceptable 
graph, it is useful to relate the latter to those of the associated metric fundamental 
polygon. The number of sides of the polygon is 2N (see above) and the corresponding 
tiling of H 2 is {p, q} with p = 2N. Due to the pairing of polygon edges, each edge of 
the graph corresponds to two edges of the polygon, i.e., p = 2e or equivalently N = e. 
In addition, each vertex of the graph corresponds to q vertices of the polygon. As the 
fundamental polygon has as many edges as vertices, p = qv. As a result, and using 
the fact that the coordinence q of the vertices in any tiling is larger than 3 (or equal), 
one finds that 2e ^ 3v, which implies that e ^ 3 (2g — 1). The only acceptable graphs 
must therefore satisfy, on top of / = 1, 

2g < e < 3 (2g - 1) 

1 < v < 2(2#- 1) 

with e and v linked by equation |6|, and g = 1 for the Euclidean plane and g ^ 2 
for the hyperbolic plane. The characteristics of the associated metric fundamental 
polygons can be inferred from the previous relations: for instance, one has 

2g < N < 3(2 5 - 1). (7) 

Note that the fact that q is an integer imposes that v must be a divider of 2e. 
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(a) (b) 



Figure 3. Constraints on the allowed closed walks on graphs: (a) each edge of 
the graph must be followed exactly once in each direction; (b) two consecutive 
steps of the walk cannot pass over the same graph edge, even if their directions 
are opposite. 



3.3. Edge pairings and closed walks on graphs 

One more step is still required. Indeed, to fully determine a fundamental polygon, the 
pairing of the edges has to be specified in addition to its shape. This pairing is closely 
related to the structure of the associated graph. We have already addressed some 
properties of the graphs: genus, numbers of vertices, edges and faces. Edge pairing 
now amounts to finding a closed walk on the graph. As explained before, going around 
the polygon boundary once is equivalent to a closed walk. This latter is constrained 
by two conditions: 

- the walk follows each edge of the graph exactly once in each direction (see 



figure 3(a)), 



- when an edge has been traversed in one direction, the walk cannot immediately 
go backward to follow the edge in the opposite direction (see figure [3(b)] ) . 

With these conditions, one can enumerate all possible graphs and closed walks for a 
given polygon. This has been done in the particular case of fundamental polygons 
with 18 sides and genus 2 [49] . 

Graphs can be represented in the Euclidean plane. However, one should remember 
that the graph does not live in this plane, but rather on the compact manifold of given 
genus g. In the following we use the term "graph" to indicate such a representation 
without a specified closed walk. With this terminology, different closed walks may 
then exist for the same graph. However, two different closed walks on a given graph 
may correspond to the same way of pairing polygon edges. For example, in the case 
of the 18-gons of genus 2 [49 , 5 graphs and 13 different closed walks were found, but 
this led to only 8 different side pairings. 

To determine the allowed closed walks on a graph, it is convenient to use the 
concept of "rotation" for vertices of the graph. For each closed walk and each vertex of 
the graph, one must specify the rules describing how one goes from one edge by which 
the walk arrives at the vertex to the subsequent one by which the walk next leaves the 
vertex, and this of course for all the edges connected to the vertex. Considering that 
for any two consecutive polygon edges their counterparts on the graph are adjacent 
edges, the only acceptable rules then correspond to clockwise and counterclockwise 
rotations around a vertex (see figure [4] for illustration). 
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Figure 4. Decorated graph and its associated clo sed w alk for the hexagonal 
periodic boundary condition in E 2 shown in figure |2(a)| As indicated by the 
arrows, the vertex with a black circle is decorated by a clockwise rotation and that 
with a white circle with a counterclockwise rotation. The present graph is a pl anar 
representation of the graph embedded in the one-hole torus (see figure [2(b)] ) . 

All possible closed walks for a given graph are then obtained after first decorating 
vertices with rotations. For each vertex there are two rotations, so that for a graph 
with v vertices there will be 2 V possible decorations. In some cases, simple rules 
allow to reduce the number of possible decorations [49] . Once all possible decorations 
have been performed, the corresponding edge pairings have to be found, knowing 
that different decorations can give the same edge pairing. A simple way to link edge 
pairings and closed walks on decorated graphs is to report on the graph the numbers 
labeling the polygon edges in the order in which they are encountered in the closed 
walk. Then, each graph edge is labeled by two numbers that indicate the pairing of 
two specific polygon sides. We now illustrate this procedure. 

Let us first detail the hexagonal case on the Euclidean plane. The graph is 
represented in figure [4j where a black vertex indicates a clockwise rotation and a 
white one a counterclockwise rotation. The displayed closed walk, which can be easily 
checked to satisfy the conditions given above, corresponds to the following edge pairing: 

1-4, 2-5, 3-6. 

This notation must be read in this way: side 1 is paired with side 4, side 2 with side 
5, side 3 with side 6. Another notation, more convenient for polygons with numerous 
sides and complex side-pairing pattern, is based on a visual description that gives, for 
each side of the polygon and when moving around in a chosen direction, the number 
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Figure 5. Decorated graph and its as sociat ed closed walk for the square periodic 
boundary condition in E 2 (see figure [l(a)] ) . The two edges form two interlaced 
closed loops; for clarity, they are labeled twice in ea ch direction. The graph 
embedded in the one-hole torus is shown in figure [T(b)1 



of sides between two paired sides. The above pairing can then be noted 
2-2-2-2-2-2, 

which means, going counterclockwise around the hexagon and starting with side 1: 1 
with 4, 2 with 5, 6 with 3. This notation will be preferred in the following. 

4. Classifying and building periodic boundary conditions 

4-1. Elements of classification 

Consider first g = 1, which corresponds to the Euclidean case. Setting g = 1 in 
equation ^ gives that the number of sides of the fundamental polygon can be 4 or 6, 
which confirms that the only two periodic boundary conditions with a regular polygon 
as primitive cell (and an associated compact quotient space) are the square and the 
hexagonal ones. We can reframe the side pairings in these two cases in the graph 
formalism of the preceding section. For the square, the graph obtained by gluing 
polygon edges and embedded in the one- hole torus has one vertex and two edges. 
So, edges form loops attached to the unique vertex. However, two arrangements are 
possible: interlaced loops and disjoint loops. By trying to perform a closed walk on 
the two possible graphs (the rotation of the unique vertex has no importance here), 
one can see that the only acceptable graph is the interlaced one shown in figure [5] 
This graph leads to the following edge pairing: 

1-1-1-1. 

For the hexagon, the graph has two vertices and three edges. Of all the 
possibilities with vertices decorated by "rotations", the only one which provides a 
suitable closed walk is that represented in figure [4] As seen previously, it leads to the 
following edge pairing: 



2-2-2-2-2-2. 
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Figure 6. The only two possible decorated graphs (and their associated closed 
walk) with one vertex and four edges. The four edges form four interlaced closed 
loops (that appear in two disconnected sets in (a), but not in (b)). Contrary to 
figure[5] the loops are only labeled once in each direction. In both cases, the genus 
is equal to 2 and the associated fundamental polygon is an octagon. 



The embedding of the graphs associated to the square and hexagonal periodic 
boundary conditions in the one-hole torus are displayed in figures [5] and [4) 

Taking now g = 2, one has to consider the hyperbolic geometry. Equation ^ for 
g = 2 becomes 

4 < TV < 9 (8) 

with, we recall, 27V denoting the number of polygon edges. For all these fundamental 
polygons, the quotient space will be a two-hole torus (see figure [9]). It is easier to 
first detail the case N = 4 for which the graph as only one vertex, like the Euclidean 
square case. Due to this unique vertex, the 4 edges become 4 loops attached to the 
vertex. In this case, only two graphs with their associated closed walk are possible 
(see figure [6]) and the corresponding side pairings are: 

1-1-5-5-1-1-5-5 and 3-3-3-3-3-3-3-3 

which correspond, respectively, to the pairing patterns shown in figures [7(a)] and [7(b)] 
In terms of the generators of the associated Fuchsian group introduced above (7^ is the 
hyperbolic translation that takes side i into its paired counterpart), the two pairings 
can also be expressed as 

7i727r 1 7^ 1 757675 _1 76 _1 = 1 and 7i7273747i" 1 7^ 1 7^ 1 7^ 1 = 1, 
where in the first case, 73 = 7^,74 = 7 2 ~ 1 ,77 = 7s" 1 ? 78 = 7^ 1 ? an d in the second 
one, 75 7f 1 >76 7 2 -1 >77 7s" 1 . 7s = 74 

Note that for g = 2 and N = 4, the fundamental polygon is an octagon of area 
equal to 47r^ -2 , which is associated to the {8,8} tiling of the hyperbolic plane. The 
link between the graph and the fundamental polygon is made by identifying paired 
polygon edges (see figure pi, which leads to the quotient space in which the graph is 
embedded (see figure [9| [36[ [43] . 

Still for g = 2, increasing N changes the shape of the fundamental polygon and its 
edge pairing, while its area remains the same. We shall not detail all the possibilities, 
and we just recall that a complete enumeration of the fundamental polygons for g = 2 
and N = 9, which corresponds to the upper bound for N in that case (see equation 
(J8|), has been given in [49] . 
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Figure 7. Fundamental polygons corr esponding to the decor ated graphs of 
figure [6] (a) corresponds to figure |6(a)| and (b) to figure |6(b)| The pairings 
are the only ones allowed for the {8, 8} tiling of H 2 . 




(a) (b) 



Figure 8. Schematic re presen tation of the "compact ificat ion" of the fundamental 
polygon shown in figure [7(a)] The paired edges are glued together: 1 with 3 and 
5 with 7 in (a); 2 with 4 and 6 with 8 in (b). The final compact manifold is a 




(a) (b) 



Figure 9. Compact manifold and graphs obtained by glu ing the edges of the 
fundamental polygons of figure [7(a)] (in(a)) and figure [7(b)] (in (b)). The graphs 
embedded in the compact manifold (a two-hole torus) correspond to the planar 
representations of figure [6] 
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Figure 10. Decorated graphs and associated closed walks of two "families" of 
periodic boundary conditions in H 2 for a generic value of the genus g ~^_2. In 
both cases, the fundamental polygon is a £>-gon with p = Ag (see figure [TT] ), and 
the associated tiling is the {4g, Ag} one. The graphs have | edges that form | 
interlaced closed loops (the loops are only labeled once in each direction). 



What happens when g increases? First, as the value of g has no upper bound, one 
finds that an infinity of periodic boundary conditions are possible in the hyperbolic 
plane. Secondly, for given g and A/", many different edge pairings are generally possible. 
There is no point in searching for an exhaustive description, and we rather focus on 
three specific "families" of edge pairings which can be easily described for any value 
of the genus g. These three families can be thought of as generalizations of the two 
Euclidean periodic boundary conditions. 

To explain the first two families, let us go back again to the Euclidean square 
case. The associated graph has only one vertex, which is related to the fact that 
2N = 4 is the minimum value for g = 1. The analogous case for g = 2, with the 
minimum value 2N = 8, has just been considered before: two different graphs, hence 
two different edge pairings, have been found. More generally, the minimum value of 
the number of polygon edges for a given genus g is 2N = 4g (see table [T]). Periodic 
boundary conditions built with such regular 4p-gons are characterized by a graph 
with a unique vertex and 2g edges, so that the vertex has a coordinence of Ag. Yet, 
multiple possibilities arise for constructing acceptable closed walks on graphs with 
such properties (only one actually for g = 1, two for g = 2, and more for higher g's). 
Two graphs, however, are simple to generate, each one with a unique possible closed 
walk on it: they are displayed in figure [To) These two graphs lead to two different side 
pairings, and this is valid for each value of g. As illustrated in figure [TT] the associated 
fundamental polygons are 4g-gons which both represent the primitive cell of the same 
underlying {4#, 4g} tiling. 

As seen previously, the connection between the "graph" and the "polygon" points 
of view is achieved by gluing polygon edges to obtain the quotient space where the 
edges now form the graph (see figure 12(a)). The two families we have constructed 
in the above described manner can be considered as generalizations of the Euclidean 
square case, which we retrieve by taking g = 1. In the case represented in figure [TT(a)| 
it is noteworthy that the metric fundamental polygon is also a standard (or canonical) 
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(a) 



(b) 



Figu re 11. Fundamental p olygons associated with the decorated graphs of 
figure |10(a)| (in (a)) and figure |10(b)| (in (b)). 




(a) 



(b) 



Figure 12. Comp act ma nifold obtained by g luing t he edges of the fundamental 
polygons of figure |ll(a)| (in (a)) and figure |ll(b)| (in (b)). The manifold is 
generically a g-hole torus; in the present case, g = 3. The graphs of figure ^0] 
are shown here embedded in the torus. 



fundamental polygon, which is not true for the family illustrated in figure |ll(b)[ 
However, from a practical point of view, there is no difference in the effort needed to 
implement these two types of periodic boundary conditions. 

Another simple "family" is made up of the hyperbolic equivalents of the Euclidean 
hexagonal periodic boundary conditions. This family corresponds to fundamental 
polygons of genus g with 27V = 4g + 2 sides (see table [I]). The main characteristic 
of this family is the relatively straightforward way to pair polygon sides. Indeed, 
in this case, pairing consists in linking opposite sides, regardless of the value of g. 
The corresponding graph has two vertices and 2g + 1 edges. Among all the possible 
graphs with such properties, the one corresponding to the "hexagonal" family is such 
that when the two vertices have different "rotations", edges do not intersect each 



others. This leads to a unique closed walk illustrated in figure 13(a) and to a unique 
side pairing for each value of g. (Note that this graph can be twisted and is then 
equivalent to the graph with identical rotations of the two vertices but with crossing 
edges displayed in figure [13(b)] the associated side pairing is however unchanged.) The 
fundamental polygon characterizing this family is a 2(2g + l)-gon which is associated 
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(b) 



Figure 13. (a) Decorated graph of a third "family" of periodic boundary 
conditions in if 2 for a generic value g ^ 2 of the genus. This family generalizes 
the hexagonal Euclidean case (g = 1) shown in figur e [4] T he corresponding 
fundamental polygon is similar to that displayed in figure [lT(b) | with p = 2(2g + l). 
The square brackets denote the integer part: here, [p/4J = g, [3p/4j = 3# + 1, 
etc. The case illustrated here corresponds to an even value of g. (For an odd 
value, one simply interchanges [p/4j and [3p/4j in the labels, so that odd labels 
always correspond to walks from left to right and opposite for even labels.) (b) 
Twisted version of the graph shown in (a), with two black vertices corresponding 
to clockwise rotations. Note that the two graphs in (a) and (b) are strictly 
equivalent in the sense that they lead to the sam e closed walk and the same 
embedded graph in the g-hole torus (see figure [l4|. 
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Figure 14. Com pact ma nifold obtained by gluing the edges of the fundamental 
polygon of figure |ll(b)| with p = 2(2g + 1), here with g = 3. The planar 
representation of the graph is shown in figure [13] 



with the underlying {2(2g + l),2g + 1} tiling and has paired opposite sides (see 



figure 11(b)). In this case, the quotient space is identical to that of the preceding 



4g-gons. It is a g-hole torus, but the embedding of the graph in the torus is different 



(see figure 14) 



4.2. Choosing the proper periodic boundary conditions 

Having sketched a classification of the possible periodic boundary conditions in the 
hyperbolic plane (see table [I]) , we now address the way to build periodic boundary 
conditions adapted to a physical problem and, therefore, satisfying some prerequisite 
constraints. For example, the symmetry of the fundamental polygon may matter in 
some cases, such as the study of ordered, crystalline-like phases on the pseudosphere. 
The surface area is also an important ingredient, e.g., when envisaging a finite-size 
analysis to study the scaling behavior of a system near a critical point. 

First, we note that due to equation Q, the area of a fundamental polygon is a 
multiple of 47r^ -2 , so that, contrary to the Euclidean case in which this area can vary 
continuously from zero to infinity, accessible areas form an infinite but countable set, 
for a given curvature. Moreover, for two different values of the area, the corresponding 
fundamental polygons must be different, since the area of a fundamental polygon 
is fixed by its genus. This sets some constraints on the associated tiling of the 
hyperbolic plane. Generically, two fundamental polygons with different areas have 
different shapes and correspond to different tilings of the hyperbolic plane. 

A fundamental polygon is characterized (in part) by its genus g and the number 
of sides 2N (recall that we only consider regular metric fundamental polygons). Each 
couple {g,N} exactly determines one tiling, {p, of the hyperbolic plane, with 



{p,q} = (2N, 2N \. 

X ^' HS \ ' N - 2g + l J 



However, different (purely hyperbolic) Fuchsian groups T of the hyperbolic plane H 2 
can produce the same tiling. This is why for the same tiling and the same shape 
of the metric fundamental polygons, the side pairings may still differ. Consider for 
example a fundamental polygon with g = 3 and 2N — 14. Different side pairings are 
possible. The simplest one is shown in figure 15(a) , whereas a more intricate one, first 
described in [47], is displayed in figure [16 (a) Yet, the two fundamental polygons are 
the primitive cell of the same {14, 7} tiling of H 2 . 

In summary, to build appropriate periodic boundary conditions on the 
pseudosphere (hyperbolic plane), one has to choose: 
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Table 1. Elements of classification of the periodic boundary conditions in the 
hyperbolic plane (conditions restricted to a regular primitive cell and a compact 
quotient space). 



Genus of 
quotient space 

Number of sides 
of the 

primitive cell 
Associated tiling 



£ii y rn.tr, 



4ff 



{45,45} 



4 5 + 2 



2N n 



Q(2g - 1) 



{2(2 5 + l),2 5 + l} {6(2 5 -l),3} 



o 



Number of edges 

Number of 
vertices 

Coordinence 
of vertices 

Number of 
possible graphs 



Allowed 
closed walks 
and decorations 



2# 
1 

> 2 C 



2^ + 1 



3(2g - 1) 



2(2g - 1) 



2^ + 1 



> 1 



May not 
exist b 

May not 
exist b 



> 5 d 





pi 
o 

b£) 
>> 

o 



pj 
S 

P 



Number of 
possible pairings 



Pairing families 



> 2 C 



> 1 



May not 
exist b 






a The Euclidean case is recovered by setting g - 
b The coordinence q of the vertices must 



1. 

q ot the vertices must be an integer and satisfy q = 
2 (1 + (2g — l)/v), where v is the number of vertices of the graph. If q is not an 
integer, the corresponding polygon and the associated tiling are not regular, which is 
a case not studied in this paper. 

c Except in the Euclidean case where the number is equal to 1. 
d The case g = 2 (see [49 ) is expected to provide a lower bound. 
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Figure 15. One possible periodic boundary condition in H 2 with g = 3 and 
2N = 14 leading to the {14,7} tiling: (a) the fundamental polygon with the 
side pairing, (b) associated compact quotient space (three-hole torus) with the 
embedded graph, and (c) decorated graph with closed walk. 



- the area of the fundamental polygon, which for a given curvature fixes its genus 
(and that of the associated quotient space), 

- the number of sides of the fundamental polygon, which for a given genus fixes 
the associated tiling (we only consider regular metric fundamental polygons and 
regular tilings), 

- the pairing of the sides, which can be achieved by using the graph formalism 
detailed above. 

From a physical point of view, the choice of a pairing among all the possible ones 
is not obvious. It might depend on relevant symmetries of the problem, but how side 
pairing affects the physics of a system is still an open question for us. On one hand, 
once periodic boundary conditions have been imposed, the boundaries of the primitive 
cell have no physical meaning. On the other hand, different side pairings correspond 
to different ways of exploring a given tiling. 




Figure 16. Same as figure [15] for another periodic boundary condition scheme 
associated to the {14, 7} tiling. Only the edge pairing patterns differ between the 
cases shown in figure [15] and figure [16] 



5. Examples of applications 

We conclude this article by illustrating how to implement the preceding developments 
and build proper periodic boundary conditions in the case of two different statistical 
mechanical models defined on the pseudosphere (hyperbolic plane). As far as we know, 
this has never been achieved before. 



5.1. The dynamics of particles on the pseudosphere 



The first example concerns the dynamics of particles on the pseudosphere. The 
starting point for implementing periodic boundary conditions is to consider the free 
motion of one particle. We pick a simple example in which the primitive cell of the 
periodic boundary condition is an octagon (genus g = 2). The trajectory of the 
particle is represented in the Poincare disk model of H 2 in figure 17 Unlike in the 



Euclidean case, this trajectory is not rectilinear. The particle under free motion follows 
geodesies, which are arcs of circles in this representation of the hyperbolic plane (see 
Appendix A). One can see in figure 7(a) the complexity of the trajectory that results 



from the pairing of the octagon sides and the resulting "folding" or "compactification" 
of space. (Alternatively, one may consider that the figure displays the trajectories of 
the particle and of all its images in the tessellation of H 2 ; the two points of view are 
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(a) 



(b) 



\ / 



i 



— --4 
(c) 



> * 

(d) 



(e) 



(f) 



(g) 



(h) 




(i) (j) (k) (1) 



Figure 17. Successive snapshots of the trajectory of a single disk-like particle 
freely evolving on the pseudosphere in the Poincare representation, with the 
periodic boundary condition shown in figure |7(a)| The trajectory follows 
geodesies, i.e., circular arcs and the arrow indicates the direction of motion. 
Note the complex pattern of exits and reentries. The fact that the particle is 
not considered as a mere point but has a finite size makes the pattern even more 
intricate. 



strictly equivalent.) A detailed account of the free motion of a mass point on the 
pseudosphere with octagonal periodic boundary conditions has been given in [36] . 

Because we intend to generalize the study to assemblies of interacting particles, 
we do not simply consider a mass point, but a particle with a given shape (here, a disk) 
and a finite size (here, a radius r*o). Note that the apparent size of the particle depends 
on the distance to the center of the Poincare disk. This stresses that the metric used 
in this representation is conformal (it preserves angles), but not isometric (it does not 
preserve lengths, which seem to contract when approaching the disk boundary): see 
|Appendi 



Let us detail a little more the practical implementation of the periodic boundary 
conditions for a such system. First, since no force acts on it, the particle follows a 
geodesic at a constant velocity. In this representation, geodesies are arcs of circles 
which meet the disk boundary at a right angle (see Appendix A), so that the 
trajectories of the particle are made of such arcs of circles. When the particle crosses 
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a side i of the octagon, the generator (a hyperbolic translation) ji linking this side to 
its paired one, 7i(i), is applied to the particle. If the particle ends up in the interior 
of the octagon, the side crossing is achieved and the particle follows another geodesic. 
But, if after this first step, the particle is still outside the octagon, one has to iterate 
the preceding step until the particle reaches the interior of the octagon. In this case, 
each new step uses the generator 7 p (i), where p is the cyclic permutation (1,2,..., 27V) 
of order 27V. Here TV = 4 as we use an octagon. Note that this procedure is valid 
only for small displacements outside the polygon, typically less than the polygon side 
length, and that the maximum number of steps is 27V — 1. 

To plot the neighborhood of the particle center, e.g., the disk of radius ro 
representing the whole particle, one duplicates the particle 27V — 1 times by using 



the preceding procedure when the disk (shown as black disk in figure 17) goes outside 
the polygon. This procedure will be necessary to study systems made of many moving 
and interacting particles. 

Implementing such periodic boundary conditions is essential to study the 
statistical mechanics of assemblies of interacting particles on the pseudosphere by 
Molecular Dynamics or Monte Carlo simulation. The physical problem under 
investigation may be crystallization and glass formation of atomic liquids in hyperbolic 
geometry p~6j [3TJ [35] . with particles modeled as disks interacting with simple pair 
potentials such as a hard-core interaction or a Lennard- Jones interaction [1 j. It 
may also be the behavior of Coulomb systems living on a pseudosphere [19] [50j [51] , 
such as the one- and two-component plasmas with interactions given by the Green's 
function associated with the Laplace-Beltrami operator. In the latter case, and as 
shown in [5j US] , the curvature introduces an "infrared cutoff' which screens the 
coulombic interactions between charged particles. As a result, the interaction is 
no longer logarithmic at large distance as it is on the Euclidean plane, but rather 
decays exponentially for distances larger than the radius of curvature ft -1 . A physical 
consequence of this curvature-induced screening is that the Kosterlitz-Thouless phase 
transition in which topological defects (vortices) all pair in dipoles is expected to be 
depressed to zero temperature. A technical consequence is that periodic boundary 
conditions can be used without recourse to specific methods such as the Ewald 
summation [2] for handling the coulombic interactions. 

In all the above cases, one is interested in studying the influence of curvature 
on the behavior of the system. The area of the primitive cell in which is embedded 
the physical system is then fixed by the genus of the fundamental polygon. However, 
there is an additional length in the problem, which is the particle size, e.g., the radius 
ro of the disk formed by an atom. Working at fixed curvature means here keeping the 
product hzro fixed. Varying the system size at constant curvature (when both the area 
and the curvature are made dimensionless by relating them to the size of one particle) 
implies to change the genus and the fundamental polygon. Keeping the dimensionless 
particle density (i.e., the "surface coverage") constant then requires to change the 
total number of particles in the primitive cell. We will not further describe how to 
implement in practice periodic boundary conditions coupled to a Molecular Dynamics 
or Monte Carlo algorithm. This will be detailed in a forthcoming publication [52] . 

5.2. Spin models on hyperbolic lattices 

We finally consider models in which spins are fixed on the vertices of a regular lattice. 
The lattices correspond to regular tilings of space, tilings which of course differ for the 
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Euclidean and the hyperbolic planes. Taking as an example of statistical mechanical 
spin system the Ising model, the Hamiltonian on a hyperbolic lattice is given by 

n = -j ^o- i0 - h 

(id) 

where &i = ±1 and the sum runs on all distinct pairs of nearest-neighbor sites, i.e., 
vertices, on the chosen lattice. 

The hyperbolic lattice Ising model has been studied in [T7\ \TE\ I20| l2T| 125] [26] 
(a continuum field theoretical version has also been considered: see [23]). We have 
already stressed in the introduction that when considering a finite-size system, surface 
effects do not become negligible compared to bulk effects when the size of the system 
increases. Indeed, the number of lattice sites around a given site grows exponentially 
with the distance, so that the contribution to any thermodynamic quantity of the 
sites at the boundary of a finite, open system is of the same order as that of the inside 
sites, even in the thermodynamic limit. Surface effects may be interesting per se [21], 
but one is nonetheless interested in "bulk" behavior, that of the "deep interior" of 
the system, because it is more directly comparable to the behavior in the absence 
of curvature. Until now, numerical simulation studies have considered finite systems 
with a free surface and studied how properties of the system evolve when neglecting 
an increasing amount of spins in the vicinity of the surface [2TJ [25] [26] . It would seem 
more convenient to study the bulk behavior by using periodic boundary conditions, 
as such conditions minimize surface effects and more directly mimic the behavior of 
an infinite system. 

Choosing the appropriate periodic boundary conditions compatible with a given 
hyperbolic lattice is however a highly nontrivial task. First, each periodic boundary 
condition corresponds to a given Fuchsian group, and so does a tiling (hence, a lattice), 
so that the Fuchsian group of the lattice and that of the periodic boundary condition 
must be compatible. Before giving a definition of what "compatible" means, it is 
useful, as before, to revisit the Euclidean case. Consider the square lattice with lattice 
spacing a. The associated discrete group of isometries has for generators the two 
translations T a $ and T a $. It is clear in this case that one can choose as a primitive cell 
any rectangle of sides mxa and nxa with m and n two nonzero integers, or restricting 
ourselves to a regular cell, any square of side nxa. The corresponding tessellation is 
also a {4,4} square tiling, and the discrete group of isometries is generated by the two 
translations T na $ and T na g. This group is a subgroup of the group associated with the 
original square lattice of spacing a: the two groups are compatible. In the hyperbolic 
case, this has to be generalized in the following way: the Fuchsian group associated 
with the fundamental polygon corresponding to the periodic boundary conditions must 
be a normal subgroup of the lattice (Fuchsian) group [44] [53] [54] . We recall that a 
subgroup T of a group E is said "normal" if it is invariant under conjugation, which 
means that for each element 7 G T and each a G E, <J7<j _1 G T. (Normal subgroups 
are closely related to quotient groups, which can indeed be constructed from a given 
group by using a normal subgroup of this group.) 

For illustrative purpose,we consider the lattice based on the {3, 7} tiling, which 
has been used for instance in previous studies of the Ising model [21] [25] [26]. The 
{3, 7} tiling is quite similar to the Euclidean triangular lattice. The normal subgroup 
compatible with the {3, 7} tiling whose fundamental polygon has the smallest genus, 
and consequently the smallest area, was first described in [47]. The metric fundamental 



polygon is a 14-gon of genus g = 3, with the side pairing displayed in figure 16(a) As 



Periodic boundary conditions on the pseudosphere 



24 




Figure 18. Periodic boundary condition adapted for a spin model on a hyperbolic 
{3, 7} lattice in the Poincare disk representation. The primitive cell in which 
the p hysical system is embedded is a 14-gon with the side pairing shown in 
figure |16(a)| Each spin on a vertex interacts with 7 nearest neighbors. The 
boundary sites shown with black circles are duplicated (see text). 



shown in figure [18] in the Poincare disk representation, it represents the unit cell of a 
{14, 7} tiling of H 2 and possesses 24 vertices of the original {3, 7} tiling. 

In figure 18, the spins which are on the boundary (black ones) are duplicated twice 
or more since they are shared between different sides of the fundamental polygon. To 
build the lattice, we have used the procedure detailed in |Appendix C For implementing 



the periodic boundary conditions, the same procedure has been followed. We have 
first constructed one cell of the {14,7} tiling associated with the normal subgroup 
of the lattice group. Then, the spins which are on the boundary of this primitive 
cell have been paired by following the fundamental polygon side pairing detailed in 



figure 16(a) Note that paired spins are only replicas of the same spin and, in the 
course of a numerical simulation, will therefore flip together (out of the 28 boundary 
spins, only 9 are independent). 

Increasing the system size, namely the size of the primitive cell, implies to choose 
another normal subgroup, and thus another fundamental polygon. The new subgroup 
corresponds to a larger value of the genus and to another tiling of H 2 (compatible 
of course with the original {3,7} tiling). The enumeration of all possible normal 
subgroups of a given Fuchsian group is still, as far as we know, an open problem [jf] 
but results have been obtained in the case of the {3, 7} tiling, and the following normal 
subgroups with larger fundamental polygons than that detailed above are described 
in [55]: their genus are equal to 7, 14, 118, 146, 411, 474, 2131 and 3404, which 
correspond, respectively, to systems with 56, 112, 944, 1168, 3288, 3792, 17048 and 
27232 spins. To further implement the corresponding periodic boundary conditions, 

|| Restrictions on the possible normal subgroups of a given Fuchsian group have been found: for 
instance in 54 , it has been shown that a {p, q} tiling is itself tiled by fundamental polygons of some 
Fuchsian group if and only if q has a prime divisor less than or equal to p. 
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one should determine the tiling associated with these subgroups and the proper pairing 
pattern for boundary spins. Work in this direction is in progress. 
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Appendix A. Hyperbolic geometry 

The hyperbolic plane H 2 , also called pseudosphere or Bolyai-Lobachevsky plane, is 
a Riemannian surface of constant Gaussian negative curvature (4TJ [42] . Contrary to 
the sphere which is a surface of constant positive curvature, H 2 as a whole cannot be 
embedded in the 3-dimensional Euclidean space E 3 . Actually, to embed H 2 in E 3 , 
the corresponding surface should have a saddle point at every point of space, which is 
obviously not feasible. The hyperbolic plane is infinite (contrary to the sphere) and 
homogeneous. 

The metric of H 2 is quite similar to that of the sphere (hence the name 
pseudosphere), particularly when expressed in polar coordinates (r, (/>); the distance 
between two infinitesimally close points is indeed given by 



where — n 2 is the negative Gaussian curvature of the plane. 

Due to the impossibility of embedding H 2 in E 3 , one has to use models to visualize 
H 2 . The representation of H 2 which may be more familiar to physicists is the upper 
sheet of an hyperboloid embedded in R 3 with a Minkowski metric. Other models 
also exist, like the Poincare upper half-plane, the Klein model, and the Poincare disk 
[41] [42] . In this work we primarily use the latter model, as it is convenient to visualize 
tilings of the hyperbolic plane. 

The Poincare disk model maps H 2 onto the open unit disk 



One can easily see that lengths in the disk are not conserved with this metric: 
the Euclidean distance between two points of A separated by a constant distance in 
H 2 shrinks to zero when the points approach the disk boundary. On the other hand, 
angles are strictly identical in H 2 and A. This representation is therefore conformal, 
but not isometric. 

The unit circle corresponding to the closure of A represents the set of points at 
infinity in H 2 . In the Poincare disk, geodesies are arcs of circles crossing orthogonally 
the limit (unit) circle. 




A = {(x,y) £ M 2 , \x 2 +y 2 \ < 1}. 
In A, the metric becomes 



_ 2 4 (dx 2 + dy 2 ) 
(l-(x 2 +y 2 )) 2 ' 



or expressed in terms of z = x + iy G C and its complex conjugate z 
4 dz dz 
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The distance p between two points z\ and Z2 is given by 



z\ - z 2 



p{zi->Z2) — 2k 1 tanh 

1 21-22 

The area of a hyperbolic disk of radius r in H 2 is equal to Attk~ 2 sinh 2 (^f ) whereas the 
length of a hyperbolic circle of radius r is given by 2ttk~ 1 sinh(^r). As a consequence, 
the ratio of the perimeter of a disk to its area goes to a nonzero value, k, when nr — > oo. 

Appendix B. Isometries and Fuchsian groups 

Physically, isometries correspond to displacements that leave the metric invariant, so 
they are an important ingredient of any study of a Riemannian space. 

As we mainly use the Poincare disk model in all this paper, we only detail the 
isometries for this representation of H 2 . The case of other representations will only 
be briefly touched on; it is well described for instance in [36] . 

The isometries of the Poincare disk can be represented as maps from C to C of 
the following form, 

az + c az + c 

cz + a ' cz + a ' 

where z,a, c G C, the overbar denotes complex conjugation, and \a\ 2 — \c\ 2 = 1 
(implying that if z is inside the unit disk A, so is its image). Note that these two maps 
correspond, respectively, to orientation-preserving and orientation-reversing linear 
fractional transformations. The group formed by the orientation-preserving maps is 
isomorphic to the quotient group SU (1, l)/{=bl}, where SU (1,1) is the pseudo-unitary 
group. These transformations can indeed be represented by complex matrices, 



of unit determinant. 

Note that for the Poincare upper half-plane, the group of orientation-preserving 
isometries is isomorphic to the projective linear group PSX(2,R) = SL(2, M)/{±1}, 
where the special linear group SL(2, R) consists of the set of 2 x 2 real matrices whose 
determinant is equal to +1. For the hyperboloid model, the group of orientation- 
preserving isometries is isomorphic to the restricted Lorentz group 6'0 + (2,l). All 
these groups are isomorphic to each other. The representations of H 2 are indeed 
all strictly equivalent, each one having its own advantages to study properties of the 
hyperbolic plane. Finally, it is useful to notice that the group of isometries is generated 
by reflections in "hyperbolic lines", i.e., geodesies. This property proves very helpful 
in practice. 

A Fuchsian group T is a discrete subgroup of the group of complex linear fractional 
transformations (also called Mobius transformations) with an invariant disc D, i.e., 
such that j(D) = D for any element 7 of T. If the unit disk A is T-invariant, then 
T can be viewed as a discrete subgroup of the group of isometries of the hyperbolic 



plane. T induces a tessellation of the plane by hyperbolic polygons (see Appendix 
|C| ). Note that Fuchsian groups preserve orientation. In the present article, we 
further restrict ourselves to "purely hyperbolic" Fuchsian groups that, except for the 
identity, do not contain transformations leaving points of invariant, such as rotations. 
(Such transformations with fixed points may be of two types, called "elliptic" and 
"parabolic"; in this context "hyperbolic" denotes a transformation that has no fixed 
points, i.e., a "translation".) 
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Appendix C. Hyperbolic trigonometry and tilings 

A regular tessellation of the hyperbolic plane by congruent, regular hyperbolic 
polygons, can be built by using the above mentioned property of the group of isometries 



(see Appendix B ), namely that this latter is generated by reflections in geodesies, which 
are inversions across circles in the Poincare disk model. For discrete subgroups, these 
geodesies correspond to the edges of the tessellation. Therefore, to construct a tiling 
of the hyperbolic plane, one just has to construct one tile and to duplicate it by using 
reflections in tile edges. 

To build an elementary polygon associated with the tiling {p, g}, where p 
denotes the number of sides of a tile and q the coordinence of the tile vertices, 
one can decompose it (since it is regular) into copies of a single triangle called the 
"orthoscheme" triangle [30] (see figure [CT] ) . For any tiling {p, g}, one has to first 
construct the appropriate triangle and then to duplicate it by using the above described 
procedure, where just one vertex of the triangle is used to generate the vertices of the 



tessellation (see figure CI). Such a triangle must have the following angles: ^, ^ and 



'- , and the duplicable vertex generating the tiling is the vertex corresponding to the 



angle (see figure CI). Constructing the triangle requires some elements of hyperbolic 
trigonometry, whose most important features are now summarized. 

First, the area A of any triangle with angles a, (3 and 7 is given by 

A = k~ 2 (n - (a + f3 + >y)) , 

where — ft 2 , we recall, denotes the Gaussian curvature of the hyperbolic plane. In 
hyperbolic geometry, the sum of the angles of a triangle is always smaller than tt and 
depends on its area, contrary to the Euclidean case where the sum of the angles is 
always equal to tt. (For the orthoscheme triangle introduced above, the sum of the 
angles is equal to n(l — ^~ 2 ^~ 2 ^~ 4 ), and is therefore less than tt due to the condition 
on hyperbolic tilings, equation ([3]).) 
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Then, for a general hyperbolic triangle with sides a, b and c and opposite angles 
a, (3 and 7, the trigonometric rules are as follows: 

sinh(fta) sinh(ftfr) sinh(^c) 

sin (a) sin(/3) sin(7) ' 

cosh(^c) = cosh(fta) cosh(nb) — sinh(^a) sinh(ftfr) cos(7), 

_ , v cos(a) cos(/3) + cos(7) 

cosh(^c) = . . , 

sin(a) sin(p) 

Note that the first two relations have Euclidean analogues, which can be recovered 
by taking k — > 0, whereas the last one has no Euclidean counterpart. This last relation 
implies that if two triangles have the same angles, then there is an isometry mapping 
one triangle onto the other. As a result, triangles are univocally determined by their 
angles. Finally, it is worth pointing out that the more familiar spherical trigonometry 
is recovered by making the replacement n — > in. 
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